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constant defined in Appendix A 
Airy function 
speed of sound 
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orifice diameter 

orifice inertial length defined by Eq. (24) 

cavity diameter 

small parameter [p (toLe) ^/Pi] ^ 

special functions defined in text 

special functions defined in text 

acoustic wave number (o)/c) 

constant defined by Eq. (7c) 

cavity depth 

resonator characteristic length 

ratio of orifice maximum velocity to speed of 
sound (q/c) 

special constants defined in Appendix A 
sound pressure amplitude 
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small parameter [Pi/p(ti)Le) defined by Eq. (32) 
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refers to difference between cavity and incident 
sound fields 

referes to resonator orifice area-averaged value 
denotes resonance 
refers to reflected wave 

denotes time differentiation 

denotes differentiation with respect to the 
independent variable, also denotes acoustic 
quantities 

denotes dimensional quantity 
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SUMMARY 


A semi- empirical fluid mechanical model of the acoustic be- 
havior of Helmholtz resonators is presented which predicts impe- 
dance as a function of the amplitude and frequency of the incident 
sound pressure field and resonator geometry. The model assumes 
that the particle velocity approaches the orifice in a spherical 
manner. The incident and cavity sound fields are connected by 
solving the governing oscillating mass and momentum conservation 
equations. The model is in agreement with the Rayleigh slug-mass 
model at low values of incident sound pressure level. At high 
values, resistance is predicted to be independent of frequency, 
proportional to the square root of the amplitude of the incident 
sound pressure field, and virtually independent of resonator geometry. 
Reactance is predicted to depend in a very complicated way upon 
resonator geometry, incident sound pressure level, and frequency. 
Nondimens ional parameters are defined that divide resonator impe- 
dance into three categories corresponding to low, moderately low, 
and intense incident sound pressure amplitudes. 

The two -microphone method was used to measure the impedance 
of a variety of resonators. The data were used to refine and verify 
the model. 
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INTRODUCTION 


Helmholtz resonators are currently being used as devices to 
control the intense sound fields generated within jet engines. Des- 
pite its extensive application, the acoustic behavior of the Helmholtz 
resonator is not well understood. As part of the NASA Lewis Research 
Center Basic Noise Research Program, a fluid mechanical model has 
been derived describing the acoustic behavior of the Helmholtz reson- 
ator as a function of resonator geometry, incident sound pressure 
level and frequency. 

Rayleigh^ credits Helmholtz as the originator of the first 
theoretical analysis describing the acoustic behavior of small cavity- 
backed resonators - appropriately called today the Helmholtz resonator. 
In a recent paper, Junger^ reviewed briefly the historical use of the 
Helmholtz resonators. Apparently, they were used by the Greeks to 
provide for a reverberation time in their open-air theaters. They 
were also used in Swedish and Danish churches as early as the thirteenth 
century. 

Rayleigh predicted the impedance of cavity-backed orifices by 
using the concept of lumped elements in a simple slug-mass mechanical 
oscillator analogy. His model is essentially non-fluid mechanical 
but predicts quite accurately resonance frequency when an empirical end 
correction is added to the slug mass. Rayleigh's derivation is quite 
straight-forward. He assumed that the mass contained in an 
orifice of thickness x* and diameter d* moves as a solid body when 
excited by an incident sound field. Crandall^ was one of the first 
to use fluid mechanical concepts to analyze the behavior of the sound 
fields near orifices. He solved for the friction resistance for the 
case of very long thin orifices (i.e., by assuming that the flow 

within the orifice is viscously fully-developed, one-dimensional, and 
excited by an oscillatory driving pressure gradient that is independent 
of the orifice axial coordinate. In a later study, Ingard"* modified 
Crandall's solution by including additional terms due to frictional 
losses over the wall containing the orifice. Based on both Crandall's 
earlier work and on further experimental and analytical studies, Ingard 
derived the following approximate prediction of the impedance of Helm- 
holtz resonators exposed to weak sound waves 
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The impedanae , as predicted by Eq. (1), is constant independent 
of the amplitude of the incident sound pressure for a specified fluid, 
sound frequency, and resonator geometry. Herein this will be called 
"linear" impedance. In 1935, Sivian^ observed that at high sound 
pressure amplitudes, the impedances of Helmholtz resonators are not 
constant but, instead, are strongly affected by the amplitude of the 



ina-ident sound pressure . Herein, this will be called "nonlinear" 
impedance. Sivian observed experimentally and confirmed theoretically 
that the acoustic resistance of resonators are, at sufficiently high 
sound pressures, linearly proportional to the acoustic particle vel- 
ocity in the orifice. Sivian* s studies prompted a variety of theoreti- 
cal and experimental studies to understand and predict the behavior 
of Helmholtz resonators at high sound pressure levels. Mellin® has 
recently reviewed the historical development of most of this work. 

In terms of predicting the nonlinear behavior of Helmholtz resonators, 
the work of Ingard and Ising^, Sirignano®, Zinn® and Hersh and 
Rogers^® are particularly important. 

Ingard and Ising conducted a detailed experimental investigation 
of the nonlinear behavior of an isolated orifice. By measuring simul- 
taneously the amplitude and phase of the sound pressure within the 
cavity and the particle velocity within the orifice (using a hot wire 
apparatus) they concluded that for low cavity sound pressure levels, 
the orifice resistance and reactance were in essential agreement with 
that predicted by Rayleigh's slug-mass model. At high cavity sound 
pressure levels, the measurements showed that the orifice resistance 
varied linearly with orifice (centerline) particle velocity; the 
measurements also showed that the orifice reactance was very insen- 
sitive to the cavity sound pressure level decreasing at the very high- 
est measured cavity sound pressure levels to roughly one-half its 
linear value. Ingard and Ising interpreted the orifice resistance 
data in terms of Bernouilli's Law suggesting that the flow behavior 
through the orifice is quasi-steady. The hot-wire measurements 
indicated that at these high sound pressure levels, the flow separates 
at the orifice forming a high velocity jet. Thus during one half 
cycle, the flow incident to the orifice is irrotational , but is highly 
rotational (in the form of jetting) after exiting from the orifice. 
During the other half of the cycle, the flow pattern is reversed. 

The loss of one-half of the reactance at these high pressure levels 
was accounted for by assuming that one half of the end correction 
is "blown" away by the exiting jet (in their experiments fVd<<l 
hence from Eq. (1) most of the reactance is due to the end correction). 
Ingard and Ising also measured the particle velocity as a function of 
axial distance from the orifice. They found that the inflow velocity 
rapidly decayed to very small values at distances of about two to 
three diameters from the orifice. 

Initially Sirignano and later Zinn, recognizing that Rayleigh's 
slug-mass model was incapable of accounting for the jetting of fluid 
from the orifice, used fluid mechanical concepts to predict the be- 
havior of the Helmholtz resonator. To simplify their models, they 
assumed that the characteristic dimensions of both the orifice and 
cavity are very much smaller than the incident acoustic wavelength 
and, further, that the acoustic flow through the orifice is one- 
dimensional, incompressible, quasi-steady, and calorically perfect. 

Both authors base their models on an integral formulation of 
the conservation of mass and momentum applied to two control volumes, 
one being the volume bounded by the orifice inlet and outlet surfaces 
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and the other the cavity. To solve these integrals, they used the 
method of successive approximations with the first order solution 
corresponding to the linear case of very small sound pressures inci- 
dent to an orifice. The orifice nonlinear behavior is introduced 
through the higher order terms and represents only a second order 
approximation to the (linear) first order solution. Thus their 
conclusions apply only to weakly nonlinear acoustic pressures and not 
to the intense sound pressures existing within rocket chambers or 
jet engines, the intended application of their models. 

Hersh and Rogers assumed that the acoustic inflow in the immediate 
vicinity of an orifice can be modelled as a locally spherical flow. 

Near the orifice, the particle flow is to a first approximation, unsteady 
and incompressible-referred to as the near field. An analytical 
model was presented of the behavior of orifices exposed to intense 
sound pressure fields which showed that orifice impedance to be 
dominated by particle flow nonlinearity. This model establishes ex- 
plicitly the quasi-steady behavior of orifices exposed to intense 
sound fields. 

The approach used by Hersh and Rogers is applied in Section 2 
below to derive a semi -empi rical fluid mechanical model of the non- 
linear acoustic behavior of Helmholtz resonators. In Section 3, ex- 
perimental data, obtained using the two-microphone method to measure 
resonator impedance, is used to both support and verify the model. 

The principal findings of the study are summarized in Section 4. 

2. FLUID MECHANICAL MODEL 

A fluid mechanical model of the behavior of Helmholtz resonators 
is described below. The model has been idealized sufficient to permit 
mathematical solution without compromising important physical character- 
istics. The derivation is divided into two categories, one correspond- 
ing to weak or moderately weak incident sound fields and the other to 
very intense incident sound fields. 

The model is based on the following key assumptions: 

(1) Sound waves approach the neighborhood of the resonator ori- 
fice in a spherical manner. The three-dimensionality of the sound 
field will account for the large changes in acoustic particle velocity 
measured by Ingard and Ising. For this assumption to be valid the 
orifice diameter mu^t be very small relative to the incident sound 
wavelength (i.e., d<<X*- orifice geometry is defined in Fig. 1). 

(2) The acoustic pressure and density are adiabatically related. 
Further, the incident sound is simple harmonic. 

(3) The amplitude of the incident acoustic pressure is small 
compared to the ambient pressure. 
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2 . 1 Approach 

Consider the special case of a sound wave approaching a resonator 
with its wavenumber k perpendicular to the plane surface containing the 
orifice (i.e., the cross-hatched surface shown in Fig. 2a). From this 
and assumption (1) above, the incident sound particle velocity field 
near the orifice is modelled as spherical and independent of the 
azimuthal angle ((). A sketch of the flow is shown entering the resona- 
tor cavity in Fig. 2b. The principal advantage of the model is its 
simplification of the governing equations of motion. Its principal 
disadvantage is its inability to accurately model the particle velocity 
entering the resonator cavity. This is illustrated in Fig. 2b wherein 
the spherical inflow model is singular at the virtual origin r=0; 
deviations from a truly spherical flow are required in order to permit 
the flow to enter the orifice in an axial manner. To avoid this sin- 
gularity, the (spherical) flow field is truncated at a hemispherical 
surface r*^=L*e defined such that at this surface the particle volume 
flow rate is equal to the actual particle volume flow rate entering 
the cavity. The consequence of this limitation is that the character- 
istic surface radius must be determined experimentally. This is 
discussed in more detail later. 

The spherical inflow model described above is assumed to be valid 
only during the half-cycles corresponding to particle inflow. The hot- 
wire measurements conducted by Ingard and Ising showed that at moderate 
and higher sound pressure amplitudes the outflow particle velocity 
field separated within the orifice forming a highly rotational jet. 

They further observed the particle flow to be symmetrical with respect 
to inflow and outflow. Thus a solution valid during the inflow half- 
cycle should via conservation of particle mass into and out of the 
cavity be valid throughout the cycle. 

The background material described above provides the basis for 
the following approach. First continuity and momentum conservation 
equations governing the motion of a harmonically driven spherically 
symmetric particle inflow are derived. Second, the resulting equations 
of motion are normalized by appropriately scaling the dependent and 
independent variables. The resulting equations are then simplified 
by retaining only the important terms. Separate equations of motion 
are then derived for the "linear" case (constant or almost constant 
impedance) and the "nonlinear" case (variable impedance). Third, 
the equations of motion are solved so as to satisfy two boundary condi- 
tions imposed by the pressure field. One is that the near field pres- 
sure must merge smoothly (asymptotically) into the harmonically oscil- 
lating incident pressure. The other is that at the hemispherical sur- 
face r=L% (see Fig. 2b), the near field pressure must be equal to the 
instantaneous cavity pressure. 

2 . 2 Derivation of Governing Equations 

The flow field is assumed to consist of a steady-state part 
containing no flow and an oscillating part. Assuming spherically 
symmetric inflow, the conservation of oscillating mass and momentum 
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may be written approximately 


0 * 1 . ' = o (2a) 

^'+ u'hH' + j_ / 1 - -!?*[— A 

di* dr* pA pv^r* ^r*/ 

-2jA' + _L__ A_ /^sin 0 AIl'') =0 

p)f2 r^2gj|^Q y 

(\-^) ^'- ^' = 0 

V p*f/ a© w a© 

where ( )' denotes fluctuating terms. 

Equation (2c) shows that the pressure gradient in the 0 
direction is proportional to y^which is quite small for most fluids 
of interest (viz. air). Hence the momentum balance in the 0 direction 
will be ignored. 

Equations (2a, b) can be simplified by proper normalization of 
the various terms and retaining only those of importance. The idea 
here is to try to anticipate the order-of -magnitude of the changes of 
the various terms in order to rank order their relative importance. 

To start with, the amplitude of the incident sound field Pi* is 
assumed to characterize the ^coustic pressure change external 
to the resonator orifice, Pi/c**’ the density change (recall that the 
acoustic pressure and density are adiabatically related) , q*the particle 
velocity change, Le the length scale and w*“'the time scale. 


(2b) 


(2c) 


The selection of the magnitude of the various normalization 
quantities above requires special comment. The hemispherical radius 
Le is chosen to characterize the scale of the displacement of the 
sound particle field near the orifice. Because of the spherical 
inflow assumption, the behavior of the particle displacement L% will 
be determined indirectly from measurements described in Section 3. 

The time scale o)*"‘is obvious as the response of the cavity is driven 
by the harmonic incident sound field. The selection of Pj* and q* is 
based on the experimental findings of Ingard and Ising. Their experi- 
mental data showed that two distinct regimes exist depending upon 
the amplitude Pj*of the incident sound field. For sufficiently low 
values of Pi* (the "linear" regime wherein orifice inertial reactance 
is much larger than orifice resistance) , 



and for sufficiently high values of i^*(the "nonlinear" regime 
wherein orifice resistance is much larger than orifice inertial 
reactance) , 
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Introducing the nondimens ional variables r, t, u, p, p 
defined as 


r-W L'e r i i u*'= i p*'. Prp i p-'= . 

0 


(4) 


the continuity and momentum equations may be written 







= O 


(5a) 
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The importance of the various terms in Eqs. (5a, b) are 
determined solely by the magnitude of their coefficients (recall that 
the dimensionless terms have been normalized to be of order unity) . 

To rank order them, the governing equations are divided into 
the linear regime where Eqs- (3a) apply and the nonlinear regime 
where Eq5. (3b) apply. 

Consider first the linear regime. The small parameters e and 6 
are introduced where 


e = 



« 





« 


(6) 


Assuming that Le* is of the order of d*, then 3<<1 for most aircraft 
applications. B is shown below to be a measure of the importance of 
fluid compressibility. Substituting these parameters into the 
continuity and momentum equations, the behavior of weak sound waves 
approaching the resonator orifice may be written 


(5b) 
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(7a) 


1E + i_ (r^u) + 4^- = O 

r r 3r 


bu 

H 


+ S u 


^u. 

3r 


(Jk)( 


VISCOUS \ _ Q 

-terms J 


(7b) 


An examination of the various terms in Eq. (7a) the continuity 
equation, shows that: 1) the second, term, which represents the 
divergence of the volume velocity, is of order unity and is by far 
the most important term; 2) the first term, which is a measure of 
the fluid compressibility, is very much smaller than the second term 
and 3) the third term, which takes into account the fluid nonlinearity, 
is the the smallest. Thus the flow field for the case of weak sound 
waves approaching the resonator orifice behaves predominately as 
unsteady and incompressible . This incompressible behavior of the flow 
field near the resonator orifice is consistent with the large changes 
in acoustic velocity measured by Ingard and Ising. Since the changes 
occurred over a distance very much smaller than the sound wavelength, 
they represent hydrodynamic rather than acoustic changes. 

For the nonlinear regime^ it is convenient to introduce the 
small parameters M and E defined as 



The restriction that M^<<1 follows from assumption (3) and Eq. (3b) 
above. The restriction that E<<1 follows from the assumption that 
the amplitude of the incident sound pressure field is large. Substi- 
tuting the parameters into the continuity and momentum equations 
yield 


M* E ■+ A-(r^u) -r M^r^u ^ = O 
df dr dr 


(9a) 


E -t- IX -I- 
dh dr 



(9b) 


Consistent with the linear regime, the continuity equation governing 
the nonlinear regime also behaves as if the sound field is unsteady 
and incompressible . Equation (9a) differs from Eq. (7a) in one interest- 
ing aspect. The third term in Eq. (9a) which is a measure of the 
importance of nonlinearity is larger by the amount 1/E than the first 
term which is a measure of the importance of compressibility. The 
converse is true for the linear case. Thus, the normalization for 
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the two regimes reveals an ordering among the various terms that is 
consistent and, therefore, lend credence to the validity of the 
approach. 

2 . 3 Boundary Conditions 


The governing equations of motion must satisfy two boundary 
conditions both of which are imposed upon the fluctuating pressure 
field during the half cycles corresponding to particle flow into 
the cavity. The boundary conditions provide the connection missing 
in previous studies between the instantaneous pressure within the 
resonator cavity and the harmonically oscillating driving pressure 
incident to the resonator orifice. One of the boundary conditions 
requires the near field pressure to merge smoothly (asymptotically) 
into the incident pressure. In dimensional terms, this boundary 
condition is written as 


Pi* ClOa) 

In nondimens ional terms, Eq. (10a) becomes 

Limi-l- p(r, -t)=cos(-t) 

P— >oo 


Using the Ingard and Ising hot-wire measurements as a guide, r^= “ 
corresponds practically to a radial distance of about three orifice 
diameters. Thus, when r*^> 3d*the oscillating pressure field will be 
assumed to be outside the near field influence of the orifice. 


The other boundary condition is considerably more complicated 
principally because of difficulties associated with the spherical 
inflow model. This boundary condition requires the fluctuating near 
field pressure at the hemispherical surface r^^= (recall that at 
this surface the particle volume flow rate equal to the 

actual volume flow rate entering the resonator cavity) to be equal to 
the instantaneous cavity pressure. The connection between the time 
pressure history of the cavity pressure and the particle flow rate 
into the cavity follows immediately from the adiabatic relationship 
between the cavity pressure and density written in dimensional terms, 

M = flla 

Si* 


where Zirlle ^q^(r=L!e >-^) is the instantaneous volume velocity flow rate 
entering the cavity volume Vt The negative sign is used because q^ 
is negative for inward directed spherical flow (see Fig. 2b). Using 
Eq. (4) to nondimensionalize Eq. (11a) yields 


SPc 

Si 


2rr L*e 

cS\f*P* 


u. (r= l,t ) 


(lib) 
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There are separate boundary conditions for the linear and ^ 
nonlinear regimes depending upon the relationships between P^^and q 
given by Eqs. (3a) and (3b). For the linear regime^ the nondimens ional- 
ized boundary condition is 


A_ pc )= - 


ZirLZc^ 


oo 




u 




( 12 ) 


For the nonlinear regime, the nondimensionalized boundary condition is 


pc ( r = I, -fc ) 


2-rrLl^c*^ 


u (r =1, -t) 


(13a) 


Equation (13a) can be rewritten in a mo re convenient form - this will 
be clear later-by replacing q'^by /Pf/p* from Eq. (3b) and using Eq . (8) 
to yield 

pc ( r =1^ t)= - E / a-n-L^eC^^^ u / r = I, t ) fl3b) 

V / 


There is one limitation to the derivation of the cavity pressure 
boundary conditions described by Eqs. (12) and (13) that must be dis- 
cussed. Recall that in their derivation, the cavity pressure responds 
both adiabatically and instantaneously (i.e., hydrodynamically) to 
the particle volume flow rate entering the cavity. This means that 
the cavity pressure cannot distinguish different spatial distributions 
of entering particle flow (e.g. spherical or axial). For application 
to the proposed spherical inflow model, angular effects in the 0 direc- 
tion associated with boundary layers over the side walls and orifice 
(see Fig. 2b) will have to be accounted for before application of 
Eqs. (12) and (13). This will be discussed in more detail in Sections 
2.4 and 2.5. 

2 . 4 The Linear Regime 

The solution to Eqs. (7a) and (7b) subject to the boundary condi- 
tions specified by Eqs. (10b) and (12 ) describe the behavior of 
Helmholtz resonators exposed to weak or moderately weak sound fields. 

In the derivation that follows, the effects of frictional losses due 
to the viscosity of the flow near and within the orifice are assumed to 
be inversely proportional to the square root of Re ( = (3Le^/v^, a char- 
acteristic length based Reynolds number). This assumption is based on 
the recent study of the acoustic behavior of orifices by Hersh and 
Rogers, which showed that the retarding action of the fluid viscosity 
acts, for a given driving excitation pressure, to decrease the magnitude 
of the acoustic velocity pumped through an orifice by an amount propor- 
tional to the local boundary layer displacement thickness along the 
orifice side walls. This is equivalent to assuming that the oscillating 
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boundary layer flow is laminar. Adoption of these assumptions will 
result in a particle inflow that is compatible with the cavity 
pressure response boundary condition defined by Eq. (12). Recall 
that in the derivation of these boundary conditions, the spherical 
particle inflow was assumed to be independent of the angle 0. 
Incorporating the above, Eq. (7b) may be rewritten 


llL i- C u +( I~£.i3^p] ^ + K u = 0 
H Jr 


Here K is a constant that will be evaluated experimentally. Equations 
(7a) and (7c) are further simplified by retaining terms to order e 
(since g<<l, the effects of compressibility would only weakly affect 
resonator behavior-hence are ignored) to yield 







(14a) 


+ S U -I- ^ + Jl_ U = 0 
H c)r Sr 


(14b) 


Integrating Eq. (14a) with respect to r yields 


u(r, -F(-fci£3Re)/r^ 


(15) 


where F(t ; £ ; Re) is an arbitrary function of integration. The negative 
sign again denotes that the particle velocity is directed inwards 
towards the resonator orifice. Substituting Eq. (15) into the 
radial momentum equation (Eq. 14b) yields 


_E_Vo 

SrKr 2r^ ^ r / 

The (*) notation denotes differentiation with respect to time. 
Integrating with respect to r yields 

-F -t- £ +- p + JL. _F = T (-t') 

^ r 


(16) 


(17) 


where T(t) is another arbitrary function of integration. From application 
of the boundary condition given by Eq. (10b), T(t)=cos(t). The boundary 
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condition given by Eq. (12) relates the instantaneous cavity pressure 
PcC^Jl^e) to the instantaneous particle velocity u(r=l,t) entering the 
orifice at r=l (recall again that to avoid the singularity at r=0 
associated with the spherical inflow assumption, the flow was arbi- 
trarily truncated at the hemispherical surface r=11 . To incorporate 
this boundary condition, Eq. (17) is differentiated with 
respect to time, r is set equal to unity (r=l) and Eq. (12) and 
T(t)=cos(t) substituted to yield 


- sin (j:^ 


(18) 


Equation (18) in its present form is misleading. It is valid 
only during the half-cycles associated with particle inflow. To 
properly model the particle velocity throughout all cycles Eq. (18) 
should be written 


F + ( . -t £ I F 1 '] F -i- F = - sin (f ) 


(18a) 


The motivation here is that during each half cycle the coordinate 
system must be switched to the appropriate side of the resonator 
orifice corresponding to particle inflow. While the. linear terms 
account for this automatically the nonlinear term.eFF does not. An 
alternate explanation is that the aoeffia-ient of F represents the 
system damping which must always be positive . 

The solution to Eq. (18a) yields directly, upon combining it 
with Eq. (15), the time and space behavior of the acoustic particle 
velocity field incident to the resonator. Because of its complexity, 
only an approximate solution will be sought. An analytic solution 
is possible, however, for the special case e = 0. 

2.4.1 Analytic Solution for e = 0 


Setting e = 0, Eq. (18a) reduces to the linear equation 

F+Ji_ F + F = - 


(19) 


Only the particular solution is of interest because the homogeneous 
solutions exhibit a time decay proportional to exp [-Kt/2/Re] , hence 
do not contribute to the steady-state behavior. 

The solution to the above equation, written in complex notation where 
it is understood that only the real part has physical meaning, is 
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( 20 ) 


F(t) = 


. it 


K 






ZTTlfe C*^' 






The particle velocity follows immediately by substituting F(t) 
above into Eq. (15) 


U (r,t 5 Re) = 


-e 


14: 




+ \ fl 


_ znL%cf^ \ 


( 21 ) 


The impedance of Helmholtz resonators is often defined in the 
literature as the (complex) ratio of the incident driving acoustic 
pressure (Pi=eit) to the "orifice area-averaged" particle velocity 
(uori) • To use this definition, it will be necessary to relate the 
particle velocity at the reference hemispherical surface of area 
(ZirLe^) to another hemispherical surface whose area is equal to the 
orifice area (Trd*V4). Using Eq. (15) which represents the conserva- 
tion of particle mass for spherical inflow, the orifice area averaged 
particle inflow is 


lA. QPI — 

2rr 

.11 

, -t ; Re) = - 8 ('JzeVe'y 


_ 2tt 


Trd 

L + J 




> 


The impedance of the Helmholtz resonator follows immediately to be 




= p(r=oo,-t:) 


Uc 


8VLV ^ 8{l>k) 


< « ^2. 
ZirL*eC 


( 22 ) 


From Eq. (22), resonance is defined to occur at zero reactance. 
This occurs at a radian frequency given by 


60 =CO^ 


res 


_ 2nL*e 


0*2 


v* 


(23) 


Equation (22) provides for a rigorous definition of the characteristic 
length Le- Setting the reactance specified in Eq. (1) equal to zero 
yields the following well known expression for the resonance frequency 
of Helmholtz resonators. 


2 _ TTd^^ 

CAJpgg 


V*de' 


(23a) 
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where 


de 


r+ 


0.85d' 


.0.&2.5(dl\ 


C24) 


is interpreted as the orifice "effective" inertia length. Combining 
Eq. (23a) and Eq. (22) and solving for the characteristic length 
leads to its definition 


= V*6o*res 
2ttc*^ 



(25) 


Thus the characteristic length Lg is independent of the cavity volume 
V*or the details of the incident sound pressure level and frequency - 
it is dependent primarily upon the orifice thickness and diameter. 

Equation (25) can be used to define an orifice discharge coeffi- 
cient (Cp) characteristic of Helmholtz resonators operating in the 
"linear" regime. The discharge coefficient is defined as the ratio 
of the orifice area-averaged particle velocity to the maximum particle 
velocity entering the cavity. This definition is equivalent to that 
used to define the discharge coefficient of orifices exposed to steady - 
state flow. Since the maximum particle inflow occurs at the hemispher- 
ical surface (r=l) the discharge coefficient is easily computed to be, 


Cr.= ^°ri = 2 7TL*e^ _ , /d*^V 

4 


(26) 


If If 

Since ds:de for most orifice geometries of practical interest to the 
aircraft industry, Cp is a small number. Mellin's review paper shows 
that low values of steady state discharge coefficients occur whenever 
the flow through the orifice is viscously controlled. This is cer- 
tainly consistent with the above derivation which shows the orifice 
resistance to be controlled by laminar boundary- layer losses. 

With Le defined by Eq. (25), the final form of the nondimensional - 
ized resonator impedance is 





) 


(27a) 


In dimensional terms, the resonator impedance is 



= 8K 





(27h) 
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where S=ird /4 is the orifice cross-sectional area. Equation (27b) 
compares favorably with Eq. (1) , the slug-mass model providing the 
constant K is appropriately defined and the radiation resistance 
term in .Eq. (1) is negligible. 

The amplitude and phase difference between the incident and 
cavity sound fields are derived by combining Eqs. (12) and (21) and 
evaluating them at r=l. Using complex notation where again only the 
real part has physical meaning and noting that 3/9t=i, then the 
instantaneous cavity pressure Pc is 



Recalling that the nondimensionalized incident sound field Px=eit 
then the amplitude and phase differences between the incident and 
cavity sound fields are 



(28a) 


and 



(28b) 


At resonance, the amplitude and phase angle differences between 
the cavity and incident sound fields simplify to 



= j(5_ and 
K 


4*ci “ decrees 


(28c) 


Thus at resonance, the amplitude difference between the cavity and 
incident sound field is due primarily to viscous boundary layer losses 
retarding the fluid motion near and within the orifice. 

To summarize, the classical Helmholtz resonator impedance derived 
originally by assuming that the particle mass within and near the ori- 
fice moves as a solid body (i.e., the lumped mass model) upon excita- 
tion by an incident sound field has been rederived from fluid mechani- 
cal conservation laws. This has accomplished two important objectives. 
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First, it establishes the validity o£ the spherical inflow model 
and second, it provides a logical model to understand the nonlinear 
acoustic behavior of the Helmholtz resonator. 

2.4.2 Weakly Nonlinear Solution 

The effects of moderate increases in incident sound pressure 
levels are modelled by Eq. ClSa^) • An approximate solution is sought 
in terms of a regular perturbation expansion in powers of e, the 
first two terms of which are 

P,(t; Re)+ ••• o (£") 


Collecting coefficients of the various powers of e, the leading term 
(e°=l) is identical to Eq. (19) while the terms proportional to e are 


F, F, ^(i%YF, = - F.lF.I 


(29) 


where loxes is defined by Eq. (23) . 


For this approach to be valid, | eFi (t ; Re) | < | F q ( t ; Re) | . This means 
that only modest deviations from the "linear" solutions described in 
Section 2.4.1 above are sought. This will place upper bound restric- 
tions on the magnitude of the amplitude of the incident sound field. 
Nonlinear solutions only for the special case aj^oj^es corresponding to 
resonance will be sought. This will simplify considerably the math- 
matical details and will provide insight regarding the effects of 
nonlinearity for the most important application of the Helmholtz 
resonator, at resonance. Substituting the real part of Eq . (20) for 
Fo(t) into the RHS of Eq. (2' 
special case of a)=o)res> 


F, (+i Re) = 


G 

■| + 4 K^1 

Re 

9 Re. 


where the (-) sign is valid : 
particle inflow and the (+) sign for the particle outflow. 


yields 

the following solution for 

the 


cos 

(30) 

9K 



L 9 rJ 


those 

half-cycles corresponding 

to 


To predict resonator impedance, the first harmonic or fundamen- 
tal component of Eq. (30) is required Using standard Fourier' analysis , 
it is straight-forward to show that the fundamental components 
are, approximately, given by 


(^j “ 2.C.OSQ:) - ^sinQfc) 

Z7 TT K 9 TT K 


(30a) 
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Combining zeroth and first-order terms, written in complex notation 
where it is understood that only the real part has physical meaning, 
yields 


e‘^ 

K 


l-£ ( ^ - I llEk) 

\Z7tt 


C31) 


For Eq. (31) to be valid, the constraint that | eF i(t) | < | F q ( t) | trans- 
lates into the following constraint on e. 


£ 


2 _ j 

27tt K 


XI ^ £ < K 

K 4 


This constraint on e can be translated directly into the following 
constraint placed on the upper bound of the amplitude of the incident 
sound field 




32 



where, from Eqs. (6) and (25) 


e=64/dl^" 






VdV 


1 ' 




(32) 


Following the same procedures used in deriving the linear 
resonator impedance defined by Eq. (27a) , the resonator orifice-area 
averaged impedance is, to order e. 


(-2) 



12 e 



(33a) 


Expressed in dimensional terms, 

(-zX.-, - 6 K ey. 


(33b) 


It is clear from Eq. (33) that the effects of moderately weak non- 
linearity are to increase only very slightly the resonator resistance 
and to decrease slightly the reactance. Following the same procedures 
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used to derive Eqs . (28a, b,c), the amplitude and phase differences 

between the cavity and incident sound fields are 


and 



C34a) 


C34bD 


From Eqs. (34a) and (34b) it is clear that both the ratio of sound 
pressure amplitude and phase angle differences between the cavity 
and incident sound fields decrease with increasing sound pressure 
Zeve Is . 

The decrease of reactance predicted by Eq. (33) is in agree- 
ment with that measured by Ingard and Ising. Ingard and Ising be- 
lieved that the decrease was caused by a reduction of effective 
orifice particle inertia length (de) related to particle flow jetting 
from the orifice. The derivation shown below supports this interpre- 
tation. Combining the reactance terms of Eqs. (27b) and (33b) yields, 
in dimensional terms, 

C-X^) - dg (35) 

V r - 9,r I 

Now assume that the effective orifice particle inertial length can 
be expanded as a power series in c so that to order e, 

(36) 


Substituting Eq. (36) for de* into Eq. (35) yields, to order e. 




^p*c*^S* = _ 4£ 




solution for (d^) , where ~ 


(d-) . _ 4L(d|') 

^ 9TT ^ 




Thus, the effective orifice inertial length (d^) is shown to decrease 
with increasing incident sound pressure level in accord with the 
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measurements o£ Ingard and Ising. The reduction of (dg) suggests 
that flow nonlinearity transfers stored inertial energy into dissi- 
pated kinetic energy - it apparently does not affect stored (cavity) 
potential energy. 

2 . 5 Nonlinear Regime 


Equations (9a) and (9b) govern the behavior of Helmholtz resona- 
tors exposed to intense sou^d fields. In this regime, the character- 
istic Reynolds number Re = qLe/v*i^ independent of frequency in con- 
trast to the Reynolds number Re = toLe^/v*characterizing the linear 
regime. Since the Reynolds number is fairly large in the nonlinear 
regime (Ingard and Ising measured particle speeds as high as 50 m/sec) i 
the contribution of the viscous terms will be ignored. 

Retaining terms only to order E, Eqs. (9a) and (9b) simplify to 


E iii 



)-o 

(37a) 

U 

dr 

II 

o 

(37b) 


The particle velocity follows immediately from Eq. (37a) to be 


u(r,tjE) = - sC-tjE) 


(38) 


where Gft;E) is an arbitrary function of integration. Substituting 
Eq. (38) into (37b) and integrating with respect to the radial coordi- 
nate r yields 


EQ -t- |G|G + p = rA) 
r 2r+ r J ^ / 


(39) 


where f(t) is yet another arbitrary function of integration. The 
absolute value is imposed upon G to insure that particle jetting 
absorbs acoustic energy throughout the cycle. 


The pressure p(r,t) must satisfy the boundary conditions speci- 
fied by Eqs. (10b) and (13b). Applying the boundary condition speci- 
fied by Eq. (lib), it is clear that £(t) = cos(t). To apply the 
boundary condition specif ied by Eq. (13b), Eq. (39) must first be 
differentiated with respect to time with Eq. (13b) substituted for 
3p/St (r=l,t) to yield. 






1 G 1<3 = — sin(-b) 


(39a) 
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Because Eq. (39a) is nonlinear, it is not possible to rOadily 
identify its resistive and reactive components. 

The solution to Eq. (39a) describes the time behavior of the 
particle velocity entering the cavity. It is not possible to seek 
approximate solutions by expanding G(t;E) in powers of E. This 
kind of approach becomes singular for values of t = (2n-l) tt/2, 
n=0, + 1, 2,... To illustrate this, the function G(t;E) is 

expanH^ed in powers of E, the first two terms of which are 

G(-fc5 e) = EG, (f)+o(E'‘) (40) 


Substituting Eq. (40) into Eq. (39a) yields the following approximate 
solutionns 



It is clear that the solutions Go(t) and Gi(t) are not valid 
near t = (2n-l)Tr/2; they are valid only near t = 2Trn. The methods 
of singular perturbation theory will be used to remove the singularity 
resulting in an approximate solution valid throughout the entire cycle 
(recall that Ingard and Ising's hot-wire measurements showed the oscil- 
lating orifice velocity to be symmetrical over each half-cycle) . The 
procedures for deriving the approximate solution are (1) the solution 
given by Eq. (39) will be assumed to be valid only over the part of 
the half-cycle defined by (4n-l/2) ir/2<t< (4n+l/2) tt/2 , (2) a second 
approximate solution will be derived valid for the part of the half- 
cycle defined by (4n-l)iT/2<t<(4n-l/2)ir/2 and (4n+l/2) Tt/2<t<(4n+l) tt/2 
and (3) a third approximate solution will be derived valid in the very 
thin (in time) region near t= (2n±l)Tr/2 . The method of singular per- 
turbation theory will be used to derive the third solution. Finally, 
with the approximate solution valid throughout the half-cycle range 
(4n- 1) T7/2^t< (4n+l) tt/2 , the Fourier coefficients ai and bi will be 
calculated. The Fourier coefficients will then be used to predict 
resonator impedance. 

Following this procedure, the details of which are summarized in 
Appendix A, the approximate solution to G(t;E) is 
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<3(-t;E) =; 


r 


]j2 COs(-k) 


4- b: 


i O+tVe) 


W*LO^^ 

>t(r) 


^ (^) (If g 

r»/4n+i)a. + t;a^ z<1L 

^ ^ 2. 4 

r.-(7)+ z‘’'<‘EV^ Y,%) ; 

„ _ C4n-f ["i T + 




2)^3 E H 


i 


C42) 


where Yo'Ch) and Yi-Cn) are special functions defined 
in Appendix A. The solution described by Eq. (42) is used to calcu- 
late the Fourier harmonic coefficients ai and bi where 

a,(E)~ 1.57-tO. I9E^'^^ + 0(E^} (43a) 


and 


b, (e)- 


0.07 - 0.43 Ln 


(e)- ^rrLeC^'^ ^^5S-2.&3 E 


(43b) 


With the harmonic components ai and bi specified, the resulting 
harmonic component of the acoustic particle velocity, written in 
complex notation where again it is understood that only the real 
part has physical meaning, is 


(tj t j E) = — 



a,(E)-i b,(E) 


C44) 


The relationship between the amplitudes and relative phase of 
the incident and cavity sound pressure fields are determined by 
combining Eqs . (44) and (13b). The ratio of the cavity-to-incident 
sound pressure amplitudes and relative phase is respectively 
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p? 



pr 


V J 


- I.57E 


27rLeC- \ 


(45a) 


and 



Recalling that E = /p^CayLp) ^/Pi it is clear that for very high inci- 
dent sound f ield s, I Pc 1 <<lPi I and (|)ci - 90 degrees. Further, the 
ratio jPc/Pil should vary almost linearly with E; the relative phase 
angle <t)ci, however, is a rather complicated function of the parameter 
E principally because of the term In (E) in the denominator. In 
evaluating Eq. (45b), the term -0.43 In (^will never become excessively 
large because of the constraint that the amplitude of the incident 
sound field be small relative to the ambient pressure. Equation (^5a) 
provides a convenient way to determine the characteristic length Le« 
For E small; 


f ^ lEl ^ ,.57 / zrrL*ec*" \ /Z 

pr ~ Pf" 


(46) 


Thus, measurements of lP(?/P^l, Pi, and of lead directly, via Eq. (46), 
to predictions of the characteristic length L g . The details of these 
measurements are presented in Section 3. For the remainder of this 
derivation, Lg will be presumed to be known. 

The resonator orifice area-averaged impedance follows immediately 
from its definition to yield 




3i -j- I b| 

8 {.L*eJ 



(47) 


where the term l/8(d/Le)^ references the particle velocity to the 
orifice inlet area ird*V4 instead of the reference hemispherical area 
2irL^^. Substitution of Eq. (43a) and (43b) for ai and bi yields 
for the resonator resistance and reactance 


(-R)orr 0-08 0.1 eH) 


(48a) 
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and 



0.05®“E[0.O7-0,43t„6-(i^J(^,S5-2.63E-.)] 


(48b) 


For E sufficiently small, Eqs . (48a) and (48b) simplify to 


a 0.08 



(49a) 


and 


O.OS|'^y E 0.01-043 /-n(E)- ^2 TTl-e.c^ ' ^ ^_^5 _ 2.63 E ^3^ 


(49b) 


Nondimensionalized with respect to the characteristic impedance pV, 
the resonator resistance and reactance are 






C50a) 


and 


ML 


a n n^f d*\ I 

(,L*ejV ^*) 


0.01 -0.43 Lr\ 


E- 


'2ElA^)f2.55-2.e3 



(50b) 


Further interpretation of the nonlinear impedance as defined 
by Eqs. (SOa) and (50b) is deferred until the experimental measurements 
leading to the definition of the characteristic length Le are des- 
cribed in Section 3. 
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3. 


EXPERIMENTAL PROGRAMS 


The two microphone method used by Dearf^and others has been 
applied to measure the impedance of the resonator geometries summarized 
in Table I. As shown in the Table, the range of the three important 
test parameters d^D* d^l!j were varied approximately one order- 

o£-magnitude . The purpose of the measurement program is two -fold. 
First, it will be used to define experimentally the parameter K 
introduced in the derivation of the linear theory (Section 2.4) and 
the characteristic length L% introduced in the derivation of the 
nonlinear theory (Section 2.5). Second, it will be used to improve 
and/or refine the model. 

The two -microphone approach and the instrumentation required is 
described in Section 3.1 below. The application of the data to the 
linear and nonlinear regimes are described in Sections 3.2 and 3.3 
respectively . 

3 . 1 Measurement System 


A schematic of the instrumentation and test set-up required to 
use the two-microphone method is shown in Figure 1. The resonator 
consists of a cylindrical cavity of diameter D* depth L* and an 
orifice of diameter d* and thickness x* 


Dean derives the following expressions for the resonator 
orifice-area averaged resistance and reactance 


and 




10 


SPL CO - S PL l dX 



X 


* 






splCO-spl^c) 
10 20 



(51a) 


(51b) 


where SPL(i) -SPL(c) represents the sound pressure level difference 
(in dB) between the incident sound field and the cavity sound field 
and <j)ci represents the corresponding phase difference. The two- 
microphone method of measuring impedance requires the simultaneous 
measurement of the incident and cavity sound pressure levels and 
phase. These measurements are obtained by flush mounting one 
microphone at the cavity base and the other flush with the wall 
containing the orifice as shown in Fig. 1. It is important to locate 
the incident microphone sufficiently far from the orifice to avoid 
near field effects (our experience shows that a separation distance 
of about 4 or 5 orifice diameters is adequate) . The microphone should 
be located sufficiently close, however, so that the separation dis- 
tance is small relative to the incident sound wavelength; this is 
necessary to insure accurate measurement of the incident sound wave 
amplitude and phase. 
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A schematic of the instrumentation used to conduct the 
experiments is shown in Figure 3. To generate incident sound 
pressure levels up to 160 dB, a JBL type 2480 driver capable of 
producing in excess of 10 watts of relatively "clean" acoustic power 
is used as the sound source. The 2" diameter driver throat is coupled 
to the test section by means of a 2" to 4" diameter exponential ex- 
pansion, JBL type H-93. Sound pressure levels in excess of 150 dB 
exceed the input capability of the GR 1560-P42 preamp. A 10 dB 
microphone Attenuator, GR Type 1962-3200 has been added, which 
extends the measurement range accordingly. 

The signal generated by the Heath lG-18 audio generator is 
amplified by the McIntosh MC2100 100 watt/channel power amplifier 
to power the JBL driver. The audio generator provides a tracking 
signal for the AD-YU Synchronous Filter and phase meter system. 

The 1036 system filters the two microphone input signals to the 
tracking signal frequency ± 2.5 Hz. The AD-YU Type 524A4 Phase Meter 
reads phase angle between the signals independent of signal amplitudes. 
The phase angle output is displayed on the AD-YU Type 2001 digital 
volt meter. A General Radio-1564 1/10 octave filter together with a 
Heath Type IM2202 DVM is used to record the output signals from each 
of the two microphones. Also the two signals are observed on a 
Tektronix 533 Oscilloscope to visually note approximate phase and 
distortion effects. 


The output of the incident microphone channel of the synchronous 
filter is used as a control voltage for an automatic level control 
amplifier. This control amplifier adjusts the drive level to the 
power amplifier in such a way as to keep the incident level constant, 
independent of frequency and amplitude response irregularities in 
the loudspeaker and tunnel. 


As a convenience, a triple ganged 5 dB per step ladder attenua- 
tor is used to simultaneously increase the power amplifier drive 
level and decrease the synchronous filter input signals so that the 
control loop of the automatic level control amplifier always has 
the same gain. This has the added advantage keeping the levels 
at the AD-YU Filter input constant for all testing levels. Since 
the AD-YU Filter displays a small amplitude-phase dependency, this 
improves accuracy as well as speed of data acquisition. A test of 
both microphones mounted flush in the wind tunnel wall showed phase 
tracking within ^ .2° over a sound pressure level range of 70-150dB. 


3 . 2 Linear Regime 


3.2.1 The e=0 Case 


At resonance, the response of the cavity to the incident sound 
pressure field was shown in Eq. C28c) to be proportional to the 
square root of the Reynolds number 


Pc^ 

- k" 


pr 

res 

1 


a = 



(28a) 
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The only unknown in the above equation is the parameter K. According- 
ly, tests were conducted to determine K experimentally. The test 
procedure was straight-forward. Only those resonator geometries that 
satisfied tjjje constraint that e<0.1 were considered Crecall that 
£ = ~ the particular value e<0.1 will be explained 

later} . The geometry of the resonators tested are summarized in 
Table II. The various values of K tabulated were determined by 
measuring |Pc/Pilres resonance for an incident sound pressure 
amplitude of 70dB. An approximate fit to the data is 


1C 2* 

1 + 3 (T'7d'5 


(52) 


Figure 4 shows |Pc/Pilres predicted to within 10% accuracy 

using Eq. (52) to predict K. 

With K specified by Eq. (52), the final approximate expression 
for resonator impedance, written in d-imensional terms, is 



Res'isfance 



* 

e 


Y^oO* 


Reactance 


(53) 


where (de/d) is defined by Eq. (24) and S = iTd^/4. Assuming radiation 
resistance to be negligible and d^D<<l, the resistance as predicted 
by the real part of Eq. (53) differs from that predicted by Eq. (1) 
by less than 10% over the range of resonator orifice geometries 
satisfying 0.5<rydCl.5. A comparison between predicted and experimen- 
tal resistance measurements is shown in Figure 5. The errors associated 
with predicting resistance from Eq. (53) arise from two separate 
sources, those associated with the parameter K and those associated 
with the ratio (di/d). As shown in Figure 5, Eq. (53) can be used to 
predict the linear value of resistance to within an accuracy of 
roughly 10%. 


Equation (53) can also be used to predict the off-resonance 
behavior of resonators. Following the procedures described in Section 
2.4.1, the ratio of the cavity- to-incident sound pressure amplitudes 
jPc/Pil 3.nd relative phase angles can be shown to be 



(54a) 


(54b) 
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where K is defined by Eq. (52) and =(?^S/'Vde* Figures 6 and 7 

summarize respectively the comparison between predictions based on 
using Eqs. (54a) and (54b) and measured data. As shown in the Figures, 
the comparison is excellent. 

3.2.2 Weakly Nonlinear Case 

Measurements of the weakly nonlinear resistance and reactance 
of five resonators are summarized in Figures 8 and 9 respectively. 

The nondimensionalized resistance data can be divided into two regimes, 
one in which the resistance is virtually independent of incident sound 
pressure amplitude (say e<0.1 - this is the reason only e<0.1 data 
was considered in determining the parameter K in Section 3.2.1) and 
the other strongly dependent upon it. For e>0.5, the resistance 
values appear to converge into a behavior dependent only upon the 
nonlinear parameter e. This is in agreement with the behavior predicted 
from the model governing equation (seeEq. 18a). The damping or loss 
term in Eq. (18a), the coefficient of F, is (K//Re+e | F | ) . At very 
low values of e, the first term dominates in accord with the data 
shown in Fig. 8. For sufficiently large values of e, the nonlinear 
term e|F| dominates and one would expect the importance of the 
individual resonator parameters (r’^jd.D*!?) to become unimportant. This 
is explicitly shown below in the following empirical curve fit to 
the data summarized in Fig. 8, 



where 8 (de/d*) ^K//Re represents the orifice area based non-dimens ional - 
ized resonator resistance at e=0. Equation (55) differs negligibly 
from the real part of Eq. (33a) , derived on the basis of a regular 
perturbation expansion for values of e<0.01. For large values of e, 
however, Eq . (55) matches considerably more accurately, the measured 

data as shown in Figure 8. 

According to Eq. (33a), the nondimensionalized reactance 
data should decrease with increasing incident sound pressure 
level from its tuned resonance value of zero at e=0 . The decrease 
should be only a function of the nonlinear parameter e. 


-X*" l - -4£ 

9rr 


(56) 


The measured reactance data is summarized in Fig. 9. Despite the 
scatter in the data which arises because of the extreme sensitivity 
of the cosine function near ninety degrees (recall that at resonance, 
4>ci =90 degrees), the nondimensionalized reactance data is roughly 
dependent only upon the nonlinear parameter e. The model 
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prediction Eq. (56) is seen in Fig. 9 to over predict the measure- 
ments. As shown in Fig. 9, a "better" fit to the data is 





(57) 


In agreement with the resistance data shown in Fig. 8, nonlinear 
effects are not important for e<0.1. 

The main findings of this study is that the onset of acoustic 
nonlinearity occurs only when the parameter e>0.1 where 
£ = 64 • Significant deviations from its linear or 

e=0 values do not occur until e>0.2. 

3 . 3 Nonlinear Regime 

The derivation of the nonlinear model described in Section 2.5 
assumed the characteristic length L® was independent of both incident 
sound pressure level and frequency. Equation (46) provides a con- 
venient way to indirectly measure its value. Rewriting it in terms 
of the resonator resonant frequency defined by Eq. (23a) , 


F? 

~ 12.56 /^‘-eY 1 

/' cOres'\ 1 


Pf 

VdV 

V OO^ J 1 

/r 


At resonance, where o)=(i)^es» the ratio |P*/Pi| becomes, only a function 
of the amplitude of the inci dent sound pressure field. Measurements 
of the ratio jPc/Pilres vs /ff(afde) ^/Pf are summarize d in Figure 10. 
Although the slopes of the |Pc/Pi|res vs /^(ifd®) ^/P± data appear to 
be constant, they do exhibit slight variations with resonator geometry. 
Comparing the data shown in Fig. 10 with Eq . (58) yields the follow- 
ing estimate of the characteristic length (L®) , 


I ^ 

t-e ~ 0.26 



(59) 


The characteristic length defined by Eq. (59) should be viewed with 
caution because of the extremely limited data base used in its 
determination. For example, the ratio d^D^<l for all of the 
resonator configurations shown in Fig. 10. 

The derivation leading to Eq. (59) assumes that the character- 
istic length (Lp/d*5 is independent of frequency. If this is true, 
then the ratio ]Pc/Pi*| should vary inversely with frequency as 
predicted by Eq. (58). To verify this, the frequency of the ^ 
resonator geometry having D=1.905 cm, L^2.54 cm, t^.051 cm, d=.091 cm 
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(the o symbol in Fig. 10) was varied from about 300 Hz to 650 Hz 
and the corresponding ratios |Pc/Pxl measured. A comparison 
between measurements of iPc/Pil and predictions based on combining 
Eqs . (58) and (59) is shown in Fig. 11. The excellent agreement 
between prediction and measurement shown verifies that (Le/d) is 
independent of both sound pressure level and frequency. With the 
characteristic length Le specified by Eq. (59), the final form of 
the resonator resistance and reactance follows immediately from 
Eqs. (50a) and (50b). Nondimensionalized with respect to ^d* they 
are 


■-R*‘ 


1.17 

1 

1 fi* .. 

0.31 


n+j-i 

9 ' 




1 + -L| 
9 ' 



r-x*i 

0.19 

0.07 - 0.43 in E-Z09f^e)f<^res\ 


- 

(2.55-2.63 


L gUvJ 

[dy[ 




(61) 


The nondimensionalized form of Eqs. (60) and (61) show explic- 
itly their dependence upon the parameter E and resonator geometry. 

A comparison between predicted and measured resistance and reactance 
is shown in Fig. 12 for five different resonator geometries. The 
reactance of each configuration is "tuned" to zero at an incident 
sound pressure level of 70 dB by proper adjustment of the sound 
frequency. The resistance, as predicted by Eq. (60), is in excellent 
agreement with the data, thus confirming the assumptions made in 
its derivation. The reactance, as predicted by Eq. (61), is also 
in close agreement with data as shown in Fig- 12. Also shown in 
Fig. 12 is the close agreement between predicted (via combining 
Eqs. (45b) and (59)) and measured relative phase between the incident 
and cavity sound fields. 

As a final verification of the model, the predicted resistance, 
reactance and relative phase variation with sound frequency is com- 
pared with measurements in Fig. 13. The incident sound pressure 
level was maintained at 140 dB and the frequency varied 
from about 300 Hz to 650 Hz. The measured phase change with fre- 
quency was larger than predicted. For example, between 300 and 650 
Hz, the measured phase increased from 76 to 97 degrees, a 21 degree 
phase change, while the corresponding predicted phase change is 
only about 9 degrees. The same trend was observed for the reactance. 
Here the reactance data varied from approximately -2 at 
300 Hz to +0.5 at 650 Hz in contrast with the 
predicted variation of -1.4 at 350 Hz to -0.06 at 650 Hz. 

Part of these differences may be related to nonsinusoidal 
wave form 6f the cavity sound field. Observations on the 
oscilloscope showed the cavity wave form to be highly distorted 
when the incident sound pressure level was 140 dB 
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(the corresponding cavity sound pressure level was 129 dB at 
404 Hz) , Oscilloscope traces of the incident sound pressure 
showed it to be reasonably sinusoidal. The measured cavity 
phase may contain contributions to the fundamental harmonic 
frequency arising from nonlinear interactions not accounted for 
in the nonlinear model. Since reactance is proportional to the 
cos (j), these nonlinear induced small errors in phase (e.g. from 
Fig. 13, +4.6% error at f = 300 Hz and -8JB% error at f = 650Hz) , 
strongly amplified near 90 degrees, may account for the differ- 
ences in predicted and measured values shown in Fig. 13. Figure 
13 also shows that both predicted (Eq. 60) and measured resis- 
tance are independent of frequency. Since the resistance is 
proportional to sin 4> , the errors in phase near 90 degrees are 
unimportant. This explains the excellent agreement between 
predicted and measured resistance. 


In reducing the data, care was taken to account for two 
serious sources of measurement error. One is phase and amplitude 
distortion of the cavity sound field due to impingement upon the 
cavity microphone by the particle velocity (see Fig. (1) during 
inflow. The second is flecture of the face plate due to the 
excessively high incident SPL's. A detailed accounting of these 
effects is summarized in Appendix B. To ensure that the relative 
phase measurements were unaffected by temperature increases due to 
the high SPL levels, a thermister probe was installed in the cavity. 
Temperature measurements showed negligible changes - less than 2 
degrees Fahrenheit for values of Pi between 70 and 160 dB. 

The experimentally determined characteristic length (Le) 
provides direct information on the value of the discharge coefficient 
Cp which connects the maximum particle velocity jetting at the 
orifice vena contracta to its orifice area-averaged value. The 
spherical inflow model forces the maximum vena contracta particle 
velocity to occur at the hemispherical surface 2irLe^. It follows, 
via the continuity of particle volume flow velocity, that 

2.nL%^ = -Trd*^ ^ Cr> = iffi = 8/^^ 

Substituting Eq. (59) for the ratio Lg/d^yields 


C 


D- 


0.545’ 



(62) 


a value remarkably close to that measured by orifices exposed to 
steady-state flow (see e.g., Fig. E-1 of the study by Rogers and Hersh^) . 

j, is clear from a comparison between Eqs. (60) and (61) that 

1R]>>|X| regardless of the frequency of the incident sound (the 
constraint that E<<1 restricts the maximum value of the frequency) . 
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Thus, from a sound absorption application viewpoint, Helmholtz resona 
tors operate as broadband absorbers in the high SPL range. To clarify 
this, consider a wall containing an array o£ Helmholtz resonators 
exposed- to intense sound. The driving sound field consists of 
both incident and reflected waves. The fraction of incident 
acoustic energy (a) absorbed by the Helmholtz resonators is a 
function of the wall impedance and may be written ^ ^ 


4 

(R'L 



(^^)ori q- 1 

2. 

(x)ori 

2 



^db 



(63) 


Equations (60J and (61) rewritten in terms of p*c* m^ be used to esti- 
mate (R)^ri/P^^n^ (X)*ori/(^’*- The quantity a=(d7lJj^ is the 
total orifice open area to cavity area ratio (for cylindrical cavi- 
ties) . It is assumed that the array of Helmholtz resonators are 
sufficiently far apart from their neighbors so that they respond 
independently. Since ni)pri/pV‘ is independent of frequency (see 
Fig. 13), and | (R) ^c*cr | >> | (X) qj. j^/ pVcr | , it is clear that the 
absorption coefficient (a) is very insensitive to frequency. It 
follows that Helmholtz resonators exposed to intense sound are 
broadband sound absorbers. This is in sharp con^ast to their be- 
havior at low sound amplitudes where generally | Fr/(?c?| << | X/(?c?| 
except near resonance. 

The use of Eqs . (60) and (61) in sound absorption calculations 
is not straight forward. The complexity arises in specifying the 
incident driving sound amplitude Pf. With regard to the above 
application of a sound wave reflecting from a wall containing an 
array of Helmholtz resonators, the proper incident driving pressure 
amplitude is Pi+Pr(Pi)* Here the reflected sound field Pr(Pi) is 
written explicitly as a function of the incident sound to emphasize 
the nonlinearity of the Helmholtz wall impedance Z'iPi + Pr) , 

Pr+Pr(Pj^ ^ Z"[pr-t Pr*(Pi!>] .64 


The quantity [Z (P j^+P;r) 1 orifice-area averaged impedance, a 

is the percent open area defined earlier and u^ and u*(Pf) are the 
incident and reflected sound particle velocities respectively (far 
from the orifice). Assuming that Pj^* and Pj(Pi) are plane waves, the 
particle velocities are related to their driving pressure fields as 

Pi' = pVu* and P;(Pi^=-pt*v<* (Pf) 
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Substituting these velocity fields into Eq. (64) and solving for 
the ratio Pr(Px)/Pi yields 


Pr(Pf) , 


P* 




■h i 

■(xXrt/ffev 


+ i 

[Win /?'!’-] 


( 65 ) 


To correctly specify pj, an iterative approach to solving Eq. (65) 
is required. For example, assume that the driving pressure is the 
incident sound pressure Pf. Substitute P| into Eq. (65) to predict 
Pr. Then replace Px by Pf+Py and substitute this value back into 
Eq. (65) to obtain a revised estimate of Pj^. Repeat this process 
until P^ becomes virtually constant. In this application, the 
correct driving sound pressure amplitude incident to the resonator 
orifice is Px'*'PrCPi)* This value of the driving pressure amplitude 
would be used in calculating the parameter E in Eqs . (60) and (61). 
With regards to interpreting the data shown in Figures 10-13, the 
proper incident driving pressure amplitude is 2P* corresponding to 
pressure doubling at the incident hard wall microphone location 
(see Fig. 1) . By locating the incident microphone far from the 
resonator orifice, the complexities described above have been avoided. 

4. CONCLUSIONS 

The acoustic behavior of Helmholtz resonators has been divided 
into three categories corresponding to (1) very weak, (2) moderately 
weak and (3) very intense incident sound pressure amplitude. The 
first two categories are characterized by a "linear" or almost 
"linear" impedance and the third by a "nonlinear" impedance. The 
"linear" and almost "linear" regimes are defined by the small 
parameter e expressed below in terms of the resonator geometry, 
incident sound pressure amplitude Pi and circular frequency as 


£ = 64 

0.85 

2 

R-^ 


\ + 0.6Z5f^ 


_ 


The quantities T*and d'^represent the orifice thickness and (circular) 
diameter respectively; D*is the (cylindrical) cavity diameter. The 
corresponding nondimensionalized orifice area- averaged resistance 
and reactance at resonance are defined empirically by Eqs. (55) and 
(57), rewritten below in a slightly more convenient form, as 
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where Cp is the Reynolds number dependent orifice discharge 
coefficient defined by Eq. (26) . Both data and model prediction show 
that the effects of weak nonlinearity are to increase resistance 
and decrease reactance. The model shows resistance to arise from 
two sources, one being viscous losses that are independent of particle 
velocity and the other convective inertial losses that are linearly 
proportional to particle velocity. The decrease in reactance arises 
from a reduction of the orifice effective inertia length (de*) . In 
accord with the conclusions of Ingard and Ising, the reduction of 
di is believed to be related to particle velocity jetting generated 
by the fluid nonlinear inertia. It is clear that, providing e 
is small, both the resistance and reactance depend weakly upon the 
sound amplitude Pi*. Significant deviations of the impedance 
from its "linear" or e=0 value do not occur until e>0.2. As an 
example, consider a cavity-backed resonator with d-1.5xl0‘^ meters 
C~0.06"), T*/d*=0.2 and (d*/D*~0.1) exposed to an incident sound field 
at 1,000 Hz. Significant deviations occur when the incident sound 
pressure level is equal to or greater than about 85 dB. 


The nonlinear regime is defined by the small parameter E 
expressed below in terms of the resonator geometry and incident 
sound field as 


E = 


0.26 



fp*(co*d*) 


X 


Using the same resonator geometry and frequency as described above, 
nonlinear effects become important when the incident sound pressure 
level is equal to or greater than about 130 dB. The corresponding 
resistance and reactance follow from Eqs. (60) and (61) rewritten 
below as 


and 



where 



■r* 


1 rw 



1.57 cp y ^ 


foc*6*\ 

r .0 


\ c* 


0 . 07 - 0.43 in E - a ( 2 . 55 -2.G3E '^3) 

(1.5 


] 


and 
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is the orifice high SPL discharge coefficient. The resistance is 
shown both theoretically and experimentally to be independent of 
frequency and proportional tp the square root of the incident sound 
amplitude. Since ] R^pV| >> 1 5(7 Helmholtz resonators absorb sound 
in a broadband manner in contrast to its narrowband sound absorption 
at low pressure levels. 
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TABLE I. SUMMARY OF RESONATOR GEOMETRIES TESTED 


D’^ 

(cm) 

(cm) 

(cm) 

(cm) 

(T^d* 

(d*/D*) 

(d*/lT) 

1.905 

1.27 

0.091 

0.025 . 

0.278 

0.048 

0.072 


1.27 

0.091 

0.051 

0.558 

0.048 

0.072 


2.54 

0.091 

0.025 

0.278 

0.048 

0.036 


2.54 

0.178 

0.102 

0.573 

0.093 

0.070 


2.54 

0.356 

0.102 

0.286 

0.187 

0.140 


2.54 

0.114 

0.051 

0.071 

0.375 

0.281 


2.54 

0.178 

0.051 

0.286 

0.093 

0.070 


5.08 

0.178 

0.051 

0.287 

0.093 

0.036 


5.08 

0.356 

0.051 

0.143 

0.187 

0.070 


3.175 

1.27 

0.091 

— 0123 

0273 

OT23 

om 


2.54 

0.178 

0.025 

0.140 

0.056 

0.070 




0.356 

0.025 

0.070 

0.112 

0.140 






0.051 

0.143 










0.102 

0.278 










0.203 

0.570 








\ 

' 

0.406 

1.140 

V 

f 
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0.714 

0.051 

0.071 

0.225 

0.281 




1 


0.102 

0.143 






\ 

t 

v] 

' 

0.203 

0.284 


f 
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f 


2.5/ 


0T7T4 

0.102 

0.143 

0.14: 

. 

0.28; 

. 


3.81 0.714 0.102 0.143 0.141 0.187 
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TABLE II. SUMMARY OF RESONATOR DATA USED TO DETERMINE PARAMETER K 



T 

L* 

D^ 

fres 

1 Pc/Pi 1 max 

K 

|Pc/Pfl 

(cm) 

(cm) 

(cm) 

(cm) 

(Hz) 

(measured) 


(predicted) 

0.178 

0.025 

1.270 

1.905 

1058 

8.97 

.552 

8.84 

t? 

0.051 

2.540 

tt 

682 

7.78 

.444 

8.02 

ri 

M 

3.810 

tt 

561 

7.12 

.440 

7.27 

n 

tt 

5.080 

tt 

476 

7.07 

.408 

6.70 

0.356 

0.025 

2.540 

tt 

1153 

17.44 

.673 

17.82 

tf 

It 

It 

3.175 

697 

13.21 

.665 

13.32 

tt 

0.051 

tt 

1.905 

1097 

16.97 

.622 

18.85 

It 

tt 

3.810 

tt 

880 

16.16 

.585 

16.88 

tr 

tt 

5.080 

tt 

746 

16.11 

.540 

15.54 

It 

tt 

2.540 

3.175 

663 

14.64 

.541 

14.13 

It 

tt 

4.445 

It 

494 

13.55 

.504 

12.20 

0.714 

0.025 

3.810 

5.080 

504 

22.96 

.691 

21.99 
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Spherical Inflow Model 
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FIGURE 2. DEFINITION OF SOUND PARTICLE COORDINATE SYSTEM 
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FIGURE .3 . TWO-MICROPHONE METHOD INSTRUMENTATION 
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FIGURE 4. COMPARISON BETWEEN PREDICTED AND MEASURED CAVITY-TO-INCIDENT SOUND 
PRESSURE AMPLITUDE RATIO AT RESONANCE IN THE LOW SPL RANGE 
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FIGURE 5. COMPARISON BETWEEN PREDICTED AND MEASURED RESONATOR 
RESISTANCE AT RESONANCE IN THE LOW SPL RANGE 
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FIGURE 6. COMPARISON BETWEEN PREDICTED AND MEASURED CAVITY-TO-INCIDENT 
SOUND PRESSURE LEVEL CHANGE IN THE LOW SPL RANGE 


42 


INC - TCAV 


I 



FIGURE 7. COMPARISON BETWEEN PREDICTED AND MEASURED 
PHASE CHANGE IN THE LOW SPL RANGE 
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FIGURE 8. COMPARISON BETWEEN PREDICTED AND MEASURED WEAKLY NONLINEAR RESONATOR RESISTANCE 
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FIGURE 13. COMPARISON BETWEEN PREDICTED AND MEASURED RESONATOR 
frequency response in the high SPL RANGE 
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APPENDIX A 


DERIVATION OF NONLINEAR FOURIER COEFFICIENTS 


A detailed derivation o£ Eq. (42) and the subsequent calcula- 
tion o£ the £undamental Fourier coe££icients ai and bi Follow. 


The derivation starts From Eq. (39a) . Recall From the dis- 
cussion Following Eq. (39a) that expansion oF the Function G(t;E) 
in a regular power series oF E leads to the approximate solutions 
Go(t) and Gi(t) deFined by Eq. (41). These solutions become 
singular at t= (2n- 1) tt/ 2 . An approximate solution to G(t;E) is 
derived below valid throughout the half-oyale (4n-l) Tr/2<t< (4n+l) 17/2 
(hereaFter n=0 without loss oF generality) . 


The derivation proceeds as Follows. The halF-cycle -7r/2<t<7r/2 
is divided into the six regions shown in Fig. A-1. Region 1 is 
deFined by the interval 0<t<ir/4, Region 2 by iT/4<t<T where T will 
be speciFied later and Region 3 by T<t<iT/2. Regions -1, -2 and -3 
are deFined by the intervals -ir/4<t<0, -T<t<iT/4, and -iT/2<t<-T 
respectively. 


It will 
the Following 


prove convenient to rescale Eq. 
transFormations , 


(39a) by introducing 


(5 (i; e) = ^ (Fj s) ^ s = i/F E 


resulting in 


2 


3 






3 



(A-l) 


(A-2) 


The absolute value operator 1 ( ) | is removed because only particle 
inFlow (g>0) is considered. Since the independent variable t is 
unaFFected by the above transFormation, the regions deFined in 
Fig. A-l still apply. An approximate solution to Eq. (A-2) is 
derived valid in Regions 1 and -1 by expanding g(t;s) as Follows, 

^(■t - 3 s) (t) + S (-t)i-- (A- 3) 


Substituting Eq. (A-3) into Eq. (A-2) and collecting coeFFicients 
oF the First two powers oF s yields 


2 



and 


2 3-3' +^3-3‘ 



CA-4) 
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Approximate solutions to E.q.(A-4) are 


3- (t) = |/cos(-t) 3 




(A- 5)' 


In deriving gi(t), the function go(t) was approximated by (l-t^/4). 
This approximation restricts the solution gi(t) to Regions 1 and -1. 

To derive approximate solutions valid in Region 2, Eq. (A-2) 
is transformed as follows, 

^(T ■yS'^ = ■, Z = JL -i 


Under these transformations, Eq. (A-2) becomes 


S 


da 



3 ^ 


(A- 7) 


Expanding g(x;s) in a regular perturbation power series of s, 

j.(-C)+Sg,(T)+o(s>) 


and substituting into Eq. (A- 7) yields the following approximate 
solutions 





(A-9) 


In the derivation of gi(T), goC't) was approximated by /T(1 -t^/ 12) . 
Here A is a constant of integration. To match the solutions solutions 
go(t)+sgi(t) at t=ir/4 to go(T)+sgi(T) at t = t7/4, the interface between 
Regions 1 and 2, the constant A ==-0.59. 


Observe that Eq. (A-9) is singular at x=0 (or t=ir/2) . It is 
important to understand that the singularity occurs because the 
power series expansion used in its derivation (see Eq. A-8) assumes 
that throughout the oyole 



(A-10) 
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It is easy to show that the perturbation expansion given by Eq. 
(A-8) generates approximate solutions that violate Eq. (A-10) . 
Consider the solution described by Eq. (A-9) . Assuming that both 
T and s are small, g(T;s) =;go(T) -/t, which upon substituting into 
Eq. (A-10) and retaining only important terms yields the following 
inequality 


S 


<< I 


Obviously, this is not satisfied for x sufficiently small. The 
failure of the regular perturbation expansion described by Eq. (A-8) 
to satisfy the inequality described by Eq. (A-10) is manifested 
by the appearance of the singularities in Eq. (A-9). The method of 
singular perturbation theory is used to derive an approximate 
solution valid in Region 3. The singularities are removed by 
appropriate rescaling of the dependent and independent variables 
so that in Region 3 the terms on the LHS of Eq. (A-10) are of the 
same order as the nonlinear term on the RHS. To accomplish this, 
the following transformations are introduced 




(A-11) 


Substituting these quantities into Eq. (A- 7) yields 


2V + S 






= — COS 






CA-12) 


Since Region 3 is near x = 0(or t=ir/2) it is reasonable to assume that 
the two highest derivatives are of equal order. This occurs when 

2.m — n = 2.K1 >-nrj=| — n 


Replacing m by 1-n in Eq. (A-12) yields, after some rearrangement. 




_^vdV 

d d Y 


„ fzirLtS^ .s 

Vto*‘ 


2.-0 2n y. 2 

- I + S • 7 


(A-13) 
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Matching the leading terms on both sides 


2.-3r\ = 0 — ^ n = , m = 1/3 


Thus the appropriate equation governing the behavior of the particle 
velocity near r=0 (Region 3 in Fig. A-1) is 


dV _2Viil+ /arLiol')V= + (s%) 


where from Eq. (A-11), 

3 (^i 0 = 


CA-15-) 


The quantity T in Fig. A-1 is identified by appropriate 
matching of the solution in Region 2 to the solution of Eq. (A-14). 
As a guide to the solution of Eq. (A-14), Eq. (A-9) is written, 
for T small, 


g(CiS) = j.(T)+S^,(r)=fr(l 


__L T^Vs _L - 0-^9 
12. / at: 



Using Eq. (A-15) to replace x by n and collecting coefficients of 
the various fractional powers of s. 






4>y 


5 

-1.18 

(2rrL%S-\^4 

1 

+ 

1 

32rf''^ 



2^ 647%J 


+ o(S^) (A-16) 


Comparing Eq. (A-16) to Eq. (A-15) 
for Y(ti;s) for proper matching, 

V(^;s)^ Vo(^)+S^3 V, {rj)^ 


Substituting Eq. (A-17) into Eq. (A-14), and collecting the coefficients 
of the first two powers of s yields 

- 2 ^ _ j (A-18) 


suggests the following expansion 
Yi (jj) (A-17) 
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CA-19) 


d!^ 


2.J_ (V, V,)=o 

dri 


Integrating with respect to n. 


dVo 

I * O 

dy 



CA-20) 


dy 


- 2 V. V, =A, 


(A-21) 


where ( )'=d/dn( ) and Aq and Ai are constants of integration. The 
solution of Eq. (A-20) can be expressed in terms of the Airy function 
Ai(p) as 

Vo (7) = 

providing the constant Aq= 0. A graph of its behavior is shown in 
Fig. A- 2. The asymptotic behavior of Yo(n) can be shown (see 
Chapter 10 of the Handbook of Mathematical Functions published by 
NBS) to be 


Limit 





5 

3Zrf/z 



(A-22) 


With Yo(n) known, Eq. (A-21) can be solved numerically to determine 
Yi(n). A plot of Yi(n) is shown in Fig. A-3. Its asymptotic be- 
havior follows immediately upon substitution of Eq. (A-22) 
into Eq. (A-21) to yield. 


Limit '^\(j}') — “A 

7^- 


IZfy 47 ^ 


-t . 


25 

64 ^ 


+ O 


(f') 


(A-23) 


The method of singular perturbation theory requires that the 
asymptotic behavior of the Region 3 function Y(n;s) merge smoothly 
into the Region 2 function g(r;s). Although mathematically, the 
matching occurs asymptotically, that is in the limit as n-*-”, numeri- 
cally, matching occurs for finite values (often less than 5) of the 
independent variable. Combining Eqs. (A-15), (A-17), (A-22) and 
(A-23) yields the following asymptotic behavior for the Region 3 
functions , 
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CA-24) 


Lirni't 






5 



W 


Z5 

64 ^ 


Comparing Eq. (A- 24) above with Eq. (A- 16) which represents the 
Region 2 solution written in terms of n, matching occurs providing 
the constant 


= 1.18 f.2rrL%c?^\ 

\ / 


CA-25) 


The functions YoCh) and Yi (n) are shown in Figures A- 2 and 
A- 3 respectively to be within one percent of their asymptotic 
behavior at n=2. This value of n is used to define Region 3. 
Substituting of n = 2 in Eq. CA-15) yields r=2s^/®. From Eq. (A-6) 
and Fig. A-1, this translates into t=T=ir/2-2s^'^ . Thus an approximate 
solution valid throughout 0 <t<tt/ 2 has been constructed. In a similar 
manner, approximate solutions can be constructed for Regions -1, 

-2, and -3. Foregoing the details, the various solutions are 
described below: 

~IL 

Region - 1 ,4 ^ t - o : 

^cos {r)-s 



Re< 


jlon — Z \ {x = z + ■b') : 




AZ 


0.59 /27TL^c*^ 


(A-27) 


Region - 3 


0^ri<Z,y= T/Ey -. 

V (y) V„ (^)+ Y, (y 


(A-28) 


The notation Y denotes the function is valid_in Region- 3; the (Y^) 
denotes Region 3. Graphs of the functions YoCb) and Yi(n) are 
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shown in Figs. (A- 2) and (A- 3) respectively. 

The above approximate solutions are valid throughout the 
half-cycle -ir/2<t<iT/2 . The corresponding fundamental harmonic 
coefficients ai and hi defined below as 


^1 (e) 
b,(E) 



TT 

Z 


(3 (^-t Cos(t)d t 


Q (i ; E) sin(t)d-t 


z 


are evaluated in the usual way. The results of the integration 
are 


a, (e) ^ 1.57 -i- 0.19 E ^3 ^ 


(A- 29) 


b, (e) E 


0.07 - 0.43 in (E)- /2ttL^c*V 2.55- 2.63 E 


(A- 3 0) 


The major contribution to ai(E) comes from the regions -•n'/4£t£'rr/4 
while the major contribution to bi(E) comes from the regions 
-iT/2<t<-Tr/4 and Tr/4<t<ir/2 . 
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FIGURE A-2. DEFINITION OF THE FUNCTIONS Yo+(n) and Yq (n) 
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APPENDIX B 


CORRECTION OF HIGH SPL DATA FOR THIN ORIFICE FACE PLATES 


At very high incident sound pressure levels (150 dB and 
higher), the measured data showed that the very thin (t=0.01") 
orifice face sheets did not correlate well with the parameter E. 

This could be attributed to any, or possibly a combination of, 
three phenomena. 

(1) Cross-modulation products of higher harmonics of the cavity 
sound field could combine with the transmitted fundamental signal, 
resulting in substantial phase shifts, even if the level of such 
cross modulation products were 20 dB or more below that of the 
fundamental. Quantitative evaluation of this effect would require 
knowledge of the amplitudes and phases of all the harmonics, as well 
as an understanding of the non-linear process which generates the 
cross modulation. Since at incident levels of 150-160 dB, where 
phase correlation is poorest, the incident field contains consider- 
able 2nd harmonic distortion, there is probably some cross modulation 
occurring in the orifice as well as within the cavity, further com- 
pounding the difficulties. 

(2) Due to the physical size of the microphones required for use 
at 70 dB incident level (1/2" diameter) it is necessary to mount the 
cavity microphone at the rear of the cavity. At low to moderate 
incident sound levels, the cavity microphone is expected to measure 
the response of the entire fluid contained within the cavity. At 
high sound pressure levels, however, the influx of fluid into the 
cavity is characterized by strong jetting. If the jet does not 
decay before reaching the microphone, a stagnation pressure velocity 
would result from the impingement of the jet on the microphone 
diaphragm, giving an altered phase angle and effective cavity sound 
pressure level. Observation of the cavity sound pressure signal 

on an oscilloscope at high incident level, reveals, for the larger 
values of d/L* an instability in the wave form on one half cycle, 
which we believe to be related to this phenomenon. 

In an attempt to quantify this effect, a small hole was 
drilled into the side of one of the cavities, into which a 1/8" 
diameter microphone was inserted. The output from this microphone 
and that from the 1/2" microphone at the back of the cavity were observ- 
ed simultaneously on the oscilloscope as the incident sound pressure 
level was increased from 90 to 155 dB. At the higher levels, the 
wave instability reappeared on the 1/2" microphone trace, but was 
absent from the 1/8" microphone trace. Although the general wave 
shape of the two microphone signals was similar and there appeared 
to be a slight advance in phase of the 1/2" microphone signal with 
increasing incident sound pressure level, the phase and amplitude 
of the 1/8" microphone signal was judged too unreliable to make con- 
clusive judgements since it was extremely sensitive to slight changes 
in position. 
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(3) It was observed that the cavities with small values of t* 
and exhibited the poorest phase correlation. This suggested 
that the face plate of the cavity was being flexed by the incident 
sound pressure, effectively modulating the volume of the cavity 
and hence increasing the cavity pressure in phase with the incident 
pressure. This would have the effect of both advancing the phase 
of the cavity sound field and increasing the cavity sound pressure 
level . 


As a test of this possiblility , the case that exhibited the 
poorest correlation, (D*= 1.25", = 0.5", d^= .036" and x^= .010) 

was remeasured with the orifice covered. At low incident levels, 
the face plate flexure was inconsequential. However, at 150 dB, 
it resulted in a phase advance of 28 degrees and a cavity pressure 
increase of 1.5 dB. 

These results were sufficiently encouraging that several of 
the other samples which had given poor phase correlations were re- 
measured using this technique. In each case the phase correlation 
was improved, but the degree of improvement depended upon the face 
plate thickness and cavity dimensions. Results of the most signifi- 
cantly improved cases are plotted in Figure B-1 for phase. 

The better improvement in phase correlation achieved with the 
thinner face plates and smaller orifice diameters suggests that 
Phenomena #2 § #3 are both affecting the phase, and to a lesser 
extent the amplitude of the cavity sound pressure. We believe that 
without properly accounting for pressure change on the face plate 
due to flow through the orifice, the plate flexing correction is 
excessive. For the case of the thin face plate, small d/D'*models, 
this excess is minimal because of low airflow due to the high 
orifice resistance. Hence the phase correlation has been improved 
for these models. 

In order to obtain a better understanding of the jet impinge- 
ment phenomenon, a resonator with the dimensions D*= .75", L* = .75", 
d*= .07", T*= .02" was tested with the standard 1/2" microphone at 
the back of the cavity and the 1/8" microphone at the side of the 
cavity. The outputs from these two microphone were fed to the AD-YU 
Phase meter and simultaneously observed on the dual trace oscilloscope. 
The measurement procedure was as follows. First, the frequency of 
the incident sound was tuned at a low level to the resonance fre- 
quency. The level was then increased until the onset of the wave 
instability (from the output from the 1/2" microphone signal) 
appeared on the oscilloscope. It was observed (see Fig. B-2) that the 
instability occurred in the trough of the oscilloscope wave at an 
incident level of 145 dB. As the level was increased, the instability 
moved forward on the wave, reaching the crest of the wave at 150 dB 
and starting down the leading edge when the system limitation for the 
particular drive frequency was reached at 153 dB. This suggests a 
sound pressure level dependent delay time for whatever is causing the 
instability at the microphone. The jet from the orifice is an 
obvious candidate, since it would travel from the orifice to the 
microphone at an average velocity related to the particle velocity 
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in the orifice. 


To support the hypothesis that the jet impingement caused the 
nonlinear phase and amplitude shifts, the average jet travel time 
from the orifice to the cavity backing was measured and correlated 
with the particle velocity u in the orifice. Since the jet enters 
the cavity at the maximum (vena contractaj value of particle velocity, 
it would originate at a time corresponding to a positive-going 
zero crossing of the cavity pressure (negative going on the scope 
trace, since the condenser microphone by its nature inverts the signal) . 
The delay of the instability by some fraction (possibly >1) of a 
wavelength on the trace could then be related to its travel time if 
we were able to identify the zero crossing associated with each 
instability. To identify this, the position in time of each 
instability was noted while the oscilloscope trace was started at 
the negative-going zero crossing. Holding the amplitude constant, 
the frequency was shifted downward 10%. As expected, the instabilities 
relative position on the wave changed; the position in time of only 
one instaiblity remained constant, as illustrated by the two sketches 
shown in Figure B-3. One may conclude that the instability which 
stayed constant in time was due to the high particle velocity which 
occurred at the beginning of the oscilloscope trace, so by knowing 
the oscilloscope sweep speed, the mean travel time of the jet from 
the orifice to the microphone at the back of the cavity may be 
determined. 

As an example, for the model checked, at an incident level of 
151 dB, the delay of the jet arrival was found, from the oscilloscope, 
to be 3/4 wavelength at 784 Hz, corresponding to a travel time of 
0.96 millisecond. Since the cavity depth is 0.75" (or 0.019 meters), 
this corresponds to a mean jet velocity of 20 meters per second. 

The peak orifice velocity, in accordance with Eq. (44), would, at 
151 dB, be 33.5 meters per second. This would indicate an average 
jet speed of 60% maximum, which is not unreasonable. If we estimate 
a linear deceleration in the cavity, this would result in a peak 
pressure level of 122.5 dB at the microphone, as opposed to a peak 
level of 147 dB for the unc.ontaminated wave, or approximately 19% 
contribution. From the appearance of the oscilloscope trace this 
seems to be a high estimate, demonstrating the approximate nature of 
the linear deceleration assumption. 

The series of sketches in Figure B-3 shows how this level 
dependent travel time can contaminate the phase data from the two 
microphone method unless precautions are taken to prevent jet im- 
pingement on the cavity microphone. Note that there is (obviously) 
one jet instability per wave length, and that these would be ex- 
pected to have the form of a somewhat unstabile positive pressure 
pulse train, whose relative position on the cavity pressure wave 
depends on the mean travel time, and hence on the incident pressure, 
the cavity depth and the decay rate of the jet. 

Since the composite signal is processed through a narrow band 
filter, only the Fourier component at the test frequency will affect 
the final results. This is displayed in Fig. B-3. It will be noted 
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that the relative phase of this signal will follow the relative 
position of the pulse train. 

Figure B-4 shows the position of the pulses at three sound 
levels for the model tested at incident levels of 145, 148, 150 
and 153 dB. Visual estimate of the amplitude of the jet pulse is 
approximately 3% of the total wave. If we estimate that the 
fundamental component is 1 % of the total, the phase of the contam- 
inated signal relative to the uncontaminated would be as shown in 
the following table, for the case illustrated in Figure B5 . 


®^^inc 

ASPLinc 

A4> 

(|)c(total) 

145 

40dB 

90“ 

.6“ 

148 

40dB 

30“ 

.3“ 

150 

40dB 

-60“ 

- .5“ 

153 

40dB 

-120“ 

- .5“ 


If the amplitude of the contaminating wave were actually 3% 
of the total, the phase angle error would be increased proportionately, 
and this fraction could be even higher for small cavity depths and 
large orifices. 

One previously unexplained measurement inconsistency can now 
be at least qualitatively understood. In some of the cavities tested 
the phase was noted to approach 90“ as the level was increased, and 
then, at very high levels, to decrease again, causing reactance from 
Dean’s model to become much too negative to be physically meaningful. 

As shown earlier, this was attributable in part to face plate flexture 
and in part to jet impingement. However, it was observed that the 
phase would occasionally go past 90“ at very high levels. This effect 
could not be explained by face plate flexing. However, we see now 
that if the level, frequency and orifice to microphone distance are 
present in the correct combination to allow the jets to arrive at 
the microphone during the negative half cycles of the cavity pressure, 
a positive phase purturbation results. 

Conclusion 


We now see that in order to obtain phase and reactance data 
with sufficient accuracy to test the theoretical model at high 
incident sound levels, some basic limiting characteristic of the 
two microphone method not previously considered in the literature 
will have to be either avoided or quantified. 

The cavity flexure problem is probably not too important, 
since the very thin unsupported orifice plates for which it appears 
to be a dominant factor are probably not feasible for practical 
application. It is also a real effect which depends on the geometry 
of the test specimen rather than the measurement procedure. 

The jet impingement and cross modulation phenomena, on the 
other hand, represent effects which are detected inside the cavity. 
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but which have no bearing on the true surface impedance of the 
cavity as seen from the incident sides. The cross modulation 
problem is unavoidable. The jet impingement problem might be 
avoided by changing the microphone location. 
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FIGURE B-1. 


COMPARISON BETWEEN PREDICTED AND MEASURED PHASE 
CHANGE ACROSS ORIFICE IN THE HIGH SPL RANGE 
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Sweep speed 


Mean Travel Time = /Sweep Speed 


FIGURE B-2. DETERMINATION OF MEAN TRAVEL TIME OF JET 
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a) Uncontaminated Wave On Oscilloscope 


^ ^ ^ 


b) Pulse Train Of Jet Arrivals; Moves To The Left 
As Incident Level Is Increased 



c) Approximate Fourier Component Of Pulse Train At Test 

Frequency; Phase Advances As Incident Level Is Increased 


FIGURE B-3. PHASE RELATIONSHIP OF FUNDAMENTAL FOURIER COMPONENT 
OF JET PULSE TO UNCONTAMINATED WAVE 
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